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Analytical description of the Day-Night neutrino asymmetry 

A. D. SUPANITSKY 1 , J. C. D'OLIVO 1 , AND G. A. MEDINA-TANCO 1 . 

1 Departamento de Fisica de Altas Energias, Instituto de Ciencias Nucleares, Universidad National 
Autonoma de Mexico, A. P. 70-543, 04510, Mexico, D. E, Mexico, 
dalesupa @ gmail. com 

Abstract: We present a new treatment of the Earth matter effects on neutrino oscillations that is valid 
for an arbitrary density profile. When applied to the study of the day-night effect on the solar neutrino 
flux it renders a simple analytical expression, which is more accurate than those derived by using the 
perturbation theory and can be extended to higher energies. 



Or 
i < 

O ■ 

to: 



> 
&\ 
(N 

q 

od 
o 
i> 
o 



Introduction 

Different types of experiments have provided com- 
pelling evidence for neutrino oscillations [6]. In 
the case of solar neutrinos the leading effects can 
be accounted by oscillations between two neutrino 
flavors, parameterized in terms of the mass square 
difference Sm 2 = m 2 — m\ and the mixing an- 
gle 6. A global fit of all the existing data gives 
Sm 2 = (7.9 - 8) x 1CT 5 eV 2 and sin 2 6 = 
0.310 — 0.315 [6], which is in good agreement with 
the results of other groups. These values belong to 
the region in the parameter space referred to as the 
Large Mixing Angle Solution (LMA). According 
to the LMA, the 8 B electron neutrinos produced in 
the Sun undergo a highly adiabatic conversion and 
are almost totally converted into the mass eigen- 
state v 2 - Then, the electron neutrino survival prob- 
ability is P(y e — > v e ) = sin 2 0. However, during 
the night solar neutrinos arriving to terrestrial de- 
tectors travel a certain distance through the Earth's 
matter, which affects the oscillations pattern. This 
leads to a partial regeneration of the electron neu- 
trino flux, a phenomenon known as the day-night 
effect. 

Matter effects on the neutrino oscillations inside 
the Earth are conveniently taken into account in 
terms of the parameter e(t) = 2EV(t)/Sm 2 , 
where V(t) = \/2GFn e {t) represents the poten- 
tial energy for v e , which comes from the charged- 
current interaction with electrons. Here, Gf is the 



Fermi constant, E is the neutrino energy, and n e (t) 
is the number density of electrons along the neu- 
trino path. In terms of the Avogadro number Na, 



e(t) 



0.019 



E 



lOMeV 



n e (t) 



Na cm 



x 10" 5 eV 



2 I 



(1) 



For the favored value of Sm 2 and the energy range 
of solar neutrinos, Earth's density is such that e -C 
1. Taking advantage of this fact, perturbation the- 
ory has been applied to derive an analytical expres- 
sion for the day-night rate asymmetry to first order 
in e [1, 5], which is valid for any density profile. 
The method simplifies the numerical calculations 
and it has been subsequently improved by means 
of a second order expansion in e [4] . In this work 
we show that a convenient alternative to the per- 
turbative approach is provided by the Magnus ex- 
pansion of the evolution operator [3] and from it 
we derive a more accurate formula for the regener- 
ation probability. 



Neutrino Oscillations in Matter 

We consider a system consisting of two neutrino 
flavors, $>f = (^e, ^n), which are related to the 
mass eigenstate, ^ m ass = (^l, ^2), according to 



(2) 
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where, 

U(0)=( C ° Sd ft (3) 
w y — sin 9 cos 9 J 

The evolution operator of the system satisfies the 
equation 

i^(t,t )=H(t)U(t,t ), (4) 

with the initial condition U(t a ,t a ) = 1. The 
hamiltonian in the mass base is given by 

-»-(!*)♦ 

Tr/,\( cos 2 9 sin 9 cos 9 \ 
V{t \ sing cos 9 sin 2 9 ) ' (5) 

and its eigenvalues are 

\±(t) = \[V(t) + 5 -^±A m (t))], (6) 

with 

A m (i) = ^ V /( e (t)-cos20)2+sin 2 20. 

(7) 

Let us now write 

U(t,t ) = V(t,t )Ur(t,t ), (8) 

/ e -ia-(t,t ) \ 

P(Mo) = e -ia+(t,to) J' 

where a±(i,i ) = Jt dt'\±(t'). The operator 
U-p(t,t a ) obeys Eq. (4) but for the Hamiltonian 
tf P (Mo) = VHt,t )[H(t) - H D (t)]V(t,t ), 
where Hn{t) = diag(\-(t), A+(t)). By expand- 
ing A T to first order in we obtain an approxi- 
mated expression for H-p with vanishing elements 
in the diagonal: 

(9) 

with <^ t = Jldt'A m (t'). 

The relevant quantity is the regeneration probabil- 
ity defined as the difference between the day and 
night probabilities, F reg (E) = P 2 ^ e {E) - sin 2 9, 

where P 2 ^ e {E) = | {v e \U{t, h)\v 2 } ?■ Here, we 
determine the evolution operator in the mass base 



from Eq. (8) by evaluating U-p in terms of the 
lowest-order Magnus approximation, U-p(t, t ) — 
cxp[— i f* Q dt'H-p(t' , to)]. Proceeding in such a 
way we get 

Fre 9 (E) = ^ sin(2J) sin 29 sin(<^ t ) 

+ sin 2 (I) cos 20, (10) 

with 

7 = sin26> J dt'V(t')cos(<pt^ t >)- (H) 

In writing Eq. (10), we assumed that the potential 
is symmetric with respect to the middle point of the 
trajectory i = (t + to)/2, which is the situation for 
a medium like the Earth, with a spherically sym- 
metric density profile. By keeping the lowest order 
terms of the expansion in I, our result for F reg (E) 
reduces to the one calculated to first order in e [4]. 

In order to make a numerical comparison of the 
different formulas, we examine the case of a neu- 
trino that crosses the Earth passing trough its cen- 
ter. For the electron density we adopt the simpli- 
fied model called mantle-core-mantle [7]. Accord- 
ing to it, n e (r) is approximated by a step function 
and the radius of the core and the thickness of the 
mantle are assumed to be half of the Earth radius. 
Accordingly, we put 

r 5.953 cm- 3 , r < R e /2 
n e (r) = N A I 

[ 2.48 cm~ 3 , i? e /2 < r < R e 

(12) 

where Rq is the radius of the Earth. 
Following Ref. [4], we introduce the function 

5 W = -F^lF^rHE) F^(E)], 

t reg{E) 

(13) 

where Freg is given by a certain (approximated) 
analytical expression, Frtg act ^ is obtained from 
the exact (numerical) solution, and 

F reg (E) = ^e(t s )sm 2 9 (14) 

is the average regeneration factor evaluated at the 
surface layer. Essentially, S represents the relative 
error of the approximated expression. 
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Figure 1 : Relative error 8 vs the neutrino energy in 
the case of a neutrino that goes through the Earth 
passing by its center, for Sm 2 = 8x 1CP 5 eV 2 and 
tan 2 6 = 0.4. The dashed line and dotted blue line 
are the first and second order approximations in e, 
respectively, and the solid red line corresponds to 
the first-order Magnus result. 

Figure 1 shows 8 as a function of the neutrino en- 
ergy for a neutrino that propagates inside the Earth 



and goes through its center. F, 



(appr) 



has been com- 



puted to first and second order in V and by means 
of the result given in Eq. (10). From the figure we 
see that the relative error for the Magnus approx- 
imation is always smaller than those correspond- 
ing to the perturbative calculations. Although it in- 
creases with energy it remains smaller than ^2% 
for the largest energies of the solar neutrinos. 

Day-Night asymmetry 

As a function of the energy the day-night asymme- 
try can be expressed as[l] 



A DN {E) 



where, 



2(cos26)F reg 

-, (15) 



1 - (cos26)(F reg - cos 26») 



(cos 26) (E) = / drf(r) 



cos 26 - e(E, r) 



{s(E,r) ~ cos26>) 2 + sin 2 26 
(16) 



Figure 2: Relative error in the Day-Night asym- 
metry as a function of the neutrino energy for a 
neutrino that propagates inside the Earth crossing 
through the center (77 = 0°). The curves corre- 
spond to the first (dashed line) and second (dot- 
ted blue line) order in e and to the Mangus re- 
sult (solid red line), for Sm 2 = 8x 10~ 5 eV 2 and 
tan 2 6 = 0.4. 



Here, f(r) is the spatial distribution function of 
the solar neutrino sources [2] and e(E,r) is de- 
termined by Eq. (1) with n e (r) now representing 
the electron density within the Sun [2]. Figures 2 
and 3 show the relative error in Adn(E) as a func- 
tion of the energy for the three approximations ex- 
amined here and a neutrino trajectory with nadir 
angle 77 = 0° (neutrino passing through the Earth 
center) and rj — 30° (neutrino passing tangent to 
the core region), respectively. We used the func- 
tion f(r) corresponding to the 8 B neutrinos and 
in both cases the smallest relative error is obtained 
with our expression for the regeneration probabil- 
ity. We also see that for all the approximations the 
relative error is smaller for r\ = 30°, which is due 
to the fact that the electron density, and therefore 
e, is smaller in the mantle region of the Earth. 

Finally, we also calculate the integrated day-night 
asymmetry, 

/>oo 

A DN = 2 dE<p v (E)(cos26)(E)F reg (E) 

JE th 



X 



! - / dEcj) v (E){cas26){E) 

1 



(F reg {E) 



- cos 26* 



(17) 
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Figure 3: Relative error in the Day-Night asym- 
metry as a function of the neutrino energy for a 
neutrino that propagates inside the Earth passing 
tangent to the core region (jj = 30°). The curves 
correspond to the first (dashed line) and second 
(dotted blue line) order in e and to the Mangus re- 
sult (solid red line), for Sm 2 — 8 x eV 2 and 
tan 2 9 = 0.4. 

where <f) v {E) is the normalized flux of solar 8 B 
neutrinos and E t h = 5 MeV is the detection en- 
ergy threshold for Super-Kamiokande and SNO. 
Figure 4 shows the relative error in Aon as a func- 
tion of the cosine of the nadir angle for the three 
approximated formulas. It can be seen that there 
are two regions: one corresponding to the propa- 
gation in the mantle, < cos?/ < a/3/2, and the 
other to the propagation in the mantle and the core, 
cos 77 > a/3/2. The relative error is practically 
constant in both regions. In the mantle it takes the 
values -1.7%, 0.07%, and -0.001% for the first 
order in e, the second order in e, and formula (10), 
respectively. In the core-mantle the corresponding 
values are -4.7%, 0.51%, and 0.13%. 

Conclusions 

In this work we have applied the Magnus expan- 
sion of the time evolution operator to find approx- 
imated analytical solutions of the system of two 
neutrino flavors coupled very weakly with matter. 
From this result we derived new expressions for the 
regeneration probability and the Day-Night asym- 
metry which give better approximations to the ex- 
act numerical results than those obtained by using 
a perturb ative approach. 
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Figure 4: Relative error in the integrated Day- 
Night asymmetry as a function of the nadir angle. 
The curves correspond to the first (dashed line) and 
second (dotted blue line) order in e and to the Man- 
gus result (solid red line), for Sm 2 = 8 x 10~ 5 eV 2 
and tan 2 9 = 0.4. 
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